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Abstract 

We studied holographic insulator /superconductor phase transition in the frame- 
work of Born-Infeld electrodynamics both numerically and analytically. First we 
numerically study the effects of the Born-Infeld electrodynamics on the phase tran- 
sition, find that the critical chemical potential is not changed by the Born-Infeld 
parameter. Then we employ the variational method for the Sturm-Liouville eigen- 
value problem to analytically study the phase transition. The analytical results 
obtained are found to be consistent with the numerical results. 



1 Introduction 



The AdS / CFT correspondence [H [21 13] which relates string theory on asymptotically AdS 
spacetime to a conformal field theory on the boundary, has been extensively applied in 
nuclear physics and condensed matter physics. In particular the related developments of 
dual gravitational models which describe insulator /superconductor phase transitions in 
condensed matter physics have been proposed in |H E]. 

The Einstein-Maxwell theory coupled to a charged scalar field would be the simplest 
model for holographic superconductor. In this model, there exists a critical temperature 
Tc [1] , below which charged AdS black holes may have non-zero scalar hair. In the context 
of AdS/CFT connection, a charged scalar field is dual to an operator which carries the 
global U{1) charge. Having non-zero scalar hair indicates that the expectation value of 
dual operator is non-zero which leads to the break of global U{1) symmetry, thus can be 
understood as a second order conductor /superconductor phase transition. The behavior 
of AC conductivity in these phases [5] make the interpretation above reliable. 

There have been many discussions concerning holographic superconductor models 
among which general relativity coupled to Maxwell field plus a charged scalar field attracts 
most attention [6] -[18] and the related Lagrangian is 

£ = -^F^ - |V^^ - iqA^^l^ - m^lV-p, F = F^.F^^'' (1) 

IIowever,the non-linear extension of the original Maxwell electrodynamics to Born-Infeld 
(BI) model [28] also caused intensive investigations [19]- [27], since BI electrodynamics 
possesses many interesting physical properties as finite total energy and the invariance 
under electromagnetic duality. In string theory, BI action can also describe gauge fields 
which arise from D-brane attached with open strings p9l |30] . The effects of BI electrody- 
namics on the holographic superconductors in the background of Ads-Schwarzschild black 
hole spacetime has been analyzed numerically in [22] and analytically in [25j.The results 
shows it is much more difficult to have scalar condensation in BI electrodynamics which 
is similar to the one of holographic superconductors in Einstein- Gauss-Bonnet gravity 
[HI HH] , where the higher curvature corrections make condensation harder. 

The authors of [lO] have constructed a model describing an insulator /superconductor 
phase transition at zero temperature, using a five dimensional AdS soliton in an Einstein- 
Maxwell- charged scalar field. The generalization of holographic insulator /superconductor 
has been carried out [311 [321 [33]. The aim of the present article is to extend holographic 
insulator/superconductor in the framework of Born-Infeld electrodynamics. The main 
purpose here is to see what effect of Born-Infeld scale parameter on the holographic 
insulator/superconductor. With both numerical and analytic method, we find a bit un- 
expected result that the effect of Born-Infeld electrodynamics is too small to be seen in 
our study, but it is consistent with [27], which state that "the critical chemical potential 
is not changed by the coupling parameter 6" . 

The paper is organized as follows. In section 2 we study holographic insulator /super- 
conductor in the framework of Born-Infeld electrodynamics in the probe limit. We numer- 
ically calculate the critical chemical potential /Zc, the condensations of the scalar operators 
{O) and charge density p. We also study the conductivity of the theory. In section 3 we 
reproduce the same result as the numerical one for the critical chemical potential /ic by 
the Sturm-Liouville analytical method. The last section is devoted to conclusions. 
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2 Numerical analysis of the critical chemical poten- 
tial /ic 

In order to construct the holographic insulator/superconductor phase transition, we con- 
sider the fixed background of five dimensional AdS soliton, e.g. in the probe limit. The 
Ads soliton metric reads 

ds"" = + r\-de + dx^ + + /(r)rfx' , (2) 

fir) = - ^ (3) 

where L is the AdS radius and tq the tip of the soliton which is a conical singularity in 
this solution. Essentially this solution can be obtained from a five dimensional AdS black 
hole solution by making use of double Wick rotations [10]. We can also set L = 1 and tq 
= 1 without the loss of generality. 

We now consider a charged complex scalar field and an electric field in the above 
background, we investigate Born-Infeld electrodynamics instead of Maxwell's therefore 
the corresponding Lagrangian density can be 

C = CBi-\V^.i^-iqA^i;\^-m^\^\^ (4) 

where is a. charged complex scalar field, Cbi is the Lagrangian density of the Born-Infeld 
electrodynamics 



^B. = ^(l-y'l + ^j. F^F.^F^'^ (5) 

The BI coupling parameter b describes the nonlinearity of electrodynamics. The above BI 
lagrangian could be reduced to Maxwell lagrangian —^F^^F^" naturally in the limit 6—7-0 
. We further assume that the BI coupling parameter h is small,since the Ads Soliton of 
background spacetimes is the solution of the Einstein equation in Einstein-Maxwell-scalar 
field theory. 

To study the insulator /superconductor phase transition we will consider the following 
ansatz HOl 



(0(r), 0,0,0), 7/> = V'(r) (6) 

The equations of motion are 



,1 , /'l , f\^' , 2&/\^,3 



and 



To seek the solutions of the equations ([7]) and ([8]), we impose the boundary condition at 
the tip r = ro and the boundary r = oo. Near the boundary, the scalar field ilj{r) and the 
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scalar potential 0(r) behave as 

^ = pe+pj' (9) 

= /i-^ (10) 

where 

A± = 2 ± V4 + m2 (11) 

is the conformal dimension of the condensation operator {O) belongs to the boundary 
field theory,and ■?/'_ and are the vacuum expectation values for the boundary operator 
{O) . fi and p are the chemical potential and the charge density,respectively. 

For simplicity we will take = — ^, which is above the Breitenlohner-Freedman 
bound El], hence ip = ^ + fy. In order to have a stable boundary theory, we impose 

the boundary condition as = or ip^ = [5j|. In this paper, we shall set ■?/'_= and 
(O) = 

At the tip, we have 

ip = ao + ai log(r - ro) + aiir - tq) + . . . , (12) 

= /3o + /3i log(r - ro) + /32(r - ro) + . . . (13) 

We can impose the Neumann boundary condition ai = /3i = to make the physical 
quantities finite. Note that the equations of motion ([71) and (IHl) have the scaling symmetry, 
{ip, 0, At, Q, b) — i- {Xip, A0, A/i, g/ A, bX^). Thus we will set g = 1 in the following discussions. 
Introducing a new variable ^ = 1/r, we rewrite the equation of motion ([7]) and ([8]) into 

and 

0" + (y + ; + (2/ - ^) - ^ (1 - fe/^V'^)'^' = (15) 

We use shooting method numerically calculate the critical potential fi^, the conden- 
sations of the scalar operators (O) and charge density p for b = 0.01. Then we change 
Born-Infeld parameter b of the theory, choosing b = 0,0.1,0.2,0.3, find that the result 
varied so small that we can neglect the effect of Born-Infeld parameter b. It is a somewhat 
surprising result to us. However, our analytic calculation support this. 

Now we holographically calculate the conductivity <j{u) by perturbing the gauge field 
|5]. To do so, we consider the following ansatz for the gauge field along x direction 

A, = A,(r)e-^'^* (16) 

Neglecting the backreaction of the perturbational field on background metric, one get the 
equation of motion of : 

b 1 2f<j>'' f'<j>'' 2f<P'r ^ Aio'Al ^^^^ 



2/ 1 + 5(_i^ _ + Z^) L r3 / r2 ' 

2Ja,A'^ 2/Ag ^ /Ag ^ 2fA'^A'' 
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Figure 1: Here we plot the behaviors of {(D) and charge density p for various values of 
chemical potential ii. Here iic = 1-888. 



Near the boundary (r — )■ oo), the behavior of the gauge field is 

...4»'.f .^i^.-.. (IS) 

the holographic conductivity is given as follows [H] 

-2i/li" , tui 

In the figure 2 we have plotted the behavior of the conductivity for different Born- 
Infield parameter b when the condensation is non-zero, we find a pole at a; — ?■ suggesting 
that we have an infinite conductivity according to the Kramers-Kronig relation as expected 
for the superconductor phase. We observed the effect of Born-Infield parameter b on the 
conductivity, the gap frequency becomes larger as Born-Infield parameter b increases. 
This is different from those of Born-Infeld electrodynamics in Schwarzschild AdS black 
hole spacetime[ 



3 Analytical calculation of the critical chemical po- 
tential fic 

From the results of the numerical calculations above, we found there exists a critical 
chemical potential Hc- When chemical potential n exceeds this critical value,the solution 
is unstable and a scalar hair will emerge which means the condensations take place. This 
is the superconductor phase. For /i < Hc,the scalar field is zero, which can be viewed as 
an insulator phase. The critical chemical potential Hc can be obtained by using Sturm- 
Liouville (SL) eigenvalue method. 

When /i < /Xc,the scalar field = 0-Eq. (fT5|) reduces to 

^' ■ ^ ■ ^2f -^)bz'(j)'A(l)' = (20) 



f z ' 2 
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Figure 2: The conductivity for operator {O) for different values of h. 



By substituting 0' to a new variable ,the equation will be transformed into a first order 
ODE of Bernoulli type, see appendix A , and the general solution is 



62 (^117^8 _ 306^4 ^ 221) _ 2bCi {77z^^ - 315z^ + 495z^ - 385) z^ 



cl 



1989 

^13 1 , ^ 

\ o -^z\ +C2 

13 9 5 32^ ' 



1155 



(21) 



where Ci and C2 are the integration constants. The boundary condition (fTOj) near the 
infinity 2; = imposes Ci = which keeps a finite value and C2 = fi- Thus, 



, , , 52 117^^- 306^4 + 221^9 

6 (z) = — ^ h u 

' 1989 ^ 



(22) 



When /i — )■ /ic ,we take {z) above into the equation (fT^ .thus the equation of motion for 
■0 becomes 



ij{z) 



b2(ll7z8_30624+22l) 
1989 



\ 



m 



z^fiz) 



\ 



+ 



- V^'(^) + = (23) 



As in [9], we introduce a trial function F {z) near the boundary z = which satisfies 

ij{z) ~ (0,)z^'F(z), (24) 
We impose the boundary condition for F{z) as -F(O) = 1 and -F'(O) = O.Thus the equation 
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of motion for F{z) is 



3/(^) _ rrA _F{± ^ ^^^^ 
4 i z^f(z) 



Following the standard procedure of Sturm-Liouville eigenvalue problem [35], we introduce 
a new function T{z) 

T{z) = (1 - z^) z^^^-^ (26) 
Multiplying T{z) to both sides of the equation of F{z) ,we have 



(Ai - 1)A 

Furthermore we define another two functions U and V as 

17 2^13 ^9\2 ^ ^^17 2;2l3 ^9^ ^2 




(27) 




+ -) -^ + -) --^+ 



17 13 9 / ^ V 17 13 9 z 



+ Ai ( ( -4 - 2^ U + ^2 1 



(28) 



= T{z)/{1 - z') (29) 
Then we can obtain the eigenvalue of /i^ by substituting U and V into the expression 

2 ^ Jo V 1 — 3Q 

In order to estimate the minimum value of /i^,we assume the trial function to be F{z) = 
1 — az^, where a is a constant. In the normal case ,the r.h.s of /i^ is an expression of a and 
other parameters already given as constants, e.g Born-Infeld parameter b in our case. But 
the problem we are dealing with is a bit different for the expression of fi^ contains /i as 
well, 

/i^ = /i^ (a, b, /i) (31) 

To avoid this difficulty ,we suppose /i on the r.h.s of the above expression to be /i at first, 
and consider /i as a new constant parameter. Then we can compute the extreme value of 
IJ,^ as we have done in the normal situation. But at this time the values of a corresponding 
to the extremal eigenvalues of yU^ are related to the parameter jl, therefore, 

/^^ = («(/^) ,b,fi) = (/i, h) (32) 
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At last ,for consistency, jl is required to be /z,and we obtain an equation for /x 

/^m = /^m(/^m>^) (33) 

where a subscript "m" is added to indicate this is an equation of extreme value fi. For 
the calculation details, see appendix B. The analytical calculation shows the minimum 
eigenvalues of fi^ and the corresponding values of a for different Born-Infeld parameter 
factor b are almost the same,which is a = 0.330 and iimin = 1-890. We see the analytical 
results fit perfectly with the numerical ones. 

4 Conclusion 

In this paper ,we considered the Born-Infeld electrodynamics in the AdS soliton back- 
ground. By numerical calculations , we found there exists a critical chemical potential 
/ic-When the chemical potential /i is lower then the critical chemical potential /Uc,the AdS 
soliton background is stable and dual field theory can be viewed as an insulator. When 
jl > He ,the background begins to be unstable and a scalar hair will be developed as the 
condensation, thus can be interpreted as the superconductor phase. So the insulator and 
superconductor phase transition appeared in the Maxwell electrodynamics also exists in 
its non-linear extension. We compute the critical chemical potential /ic both analytically 
and numerically and the results of these two methods fit perfectly good. We also found 
the effect of different Born-Infeld parameters are very small that even can be neglected. 
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A Bernoulli type first order ODE 

The standard Bernoulli type ordinary equation takes the form 

y' = fit, y{t)) = g{t)y{t) + h{t)y{tr n ^ 1 (34) 
We make the substitution ,^ = 0' in the equation (1201) which reduces to 



+ -\i-(2f-'-f]hz'e (35) 



dz \f zj V 2 

It is a Bernoulli type ordinary equation and the idea to solve the equation is to 
convert it into a linear ODE. Make the substitution 

vit) = y{tf-- (36) 

thus 

v{t)' = {l-n)y{t)-'^y{t)' (37) 
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yity 




v{ty 



(38) 



so 



v{t)' = (1 - n)g{t)v{t) + (1 - n)h{t) 



(39) 



which is a standard hnear ODE, and can be solved by the standard procedure. 

B Computation of jimin 

Substituting U y into expression (130|) ,we get 



A{a) = -25m2b^Xi [X]^ + 141Af + 8953Af + 337449Ap + 8370439Af + 
142959663A,^ + 1704340979A^ + 13997457507AJ + 75168485084A^+ 
226063117896Af + 135629702768A^ - 1262172448656Af- 
2753135394624A,^ + 1021968576000A, + 2540624486400) (a^ (64A,^+ 
2848Af + 53040A^ + 517088A^ + 2740052A,2 + 7415574A, + 7977879) - 
2a (64A,^ + 2912A,^ + 55888Af + 570128Af + 3209404A,2 + 9320246Ai+ 
10774995) + 64A,^ + 2976A,^ + 58864Af + 628992Af + 3790660A,2+ 
12084582Ai + 15758847) 

C{a) = 3956121(A, - 2)(A, - 1)A,(A, + 1)(A, + 8)(A, + 9)(A, + 10)(A, + 11)(A, + 12) 
(A, + 13) (Ai + 14) (Ai + 15) (A, + 16) (A^, + 17) (A^ + 18)(2Ai + 7)(2A, + 9) 

(2Ai + 11)(2A, + 13)(2A, + 15)(2A,, + 17)(2A,, + 19) ( - ^ + 



A{a)fi + B{a) 



Xi + 1 Xi \-l 
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B{a) = {SX^ + ISOA^ + 1318Ai + 3135) [WX]^ + 1504Af + 61872A,"+ 
1461236Af + 21801935A^ + 212273937A,^ + 1333101101AJ+ 
4996108013A^ + 8010950487A,^ - 11997217421Af - 68501188673A,^- 
66321743349A,2 + 59135273262Ai + 73109116080) (2 (-8 (646^- 
35605089) Af - 96 (3686^ - 85715955) A^ - 16 (598726^- 
7536410505) A^* - 48 (2584726^ - 19251803493) A,^- 
1728 (392696^ - 1933664031) A, + 3956121A,^ + 3828259169280) + 
3956121 (Af + 70Xt + IMOX^ + 26600A,2 + 180384Ai + 483840) m^) - 
2aAi (16A," + 1632Af + 73016A,^ + 1880616A,^ + 30716889AJ+ 
329379372A,^ + 2305033064Af + 9928719744Ai^ + 21299033680A,^- 
4869591744A,2 - 119590571520A, - 169885900800) (3956121 (8A^+ 
676AJ + 24482 A^ + 495859 A,^ + 6136697A^ + 47461414A,^+ 
223683648A,2 + 586275651Ai + 652759965) (2A,(Ai + 2) + (A^ + l)m^) - 
1286^ (64A^ + 5408A^ + 195792AJ + 3929640Af + 46829160Af+ 
322610184Af + 1147216028Af + 1425431605A2 - 1194241044Ai- 
1751976837)) + A^ (l6Af + 1792A,^^ + 88864Af + 2567628A^+ 
47829603Af + 599455521AJ + 5108576957A,^ + 29104930589A,^+ 
104893992911A^ + 206701115187Af + 113193419679A2- 
236408070717A, - 223243908030) (3956121 (8A,^ + 612AJ + 19974A,^+ 
362715Af + 4001532A^ + 27406788A,^ + 113514256A,2+ 
259082880Ai + 248371200) (2Ai + m^) - 1286^ (64A^ + 4832AJ+ 
154832A^ + 2705144A,^ + 27068168Af + 140926760Af+ 
253029340A,2 - 493864887Ai - 1682001090)) 

Taking /i = /i, Aj = | , b=0.01 and solving the above equation, we get 

flmin = 1.89037 

where we have discarded the negative value of fimin- 
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